Friction between Polymer Brushes 
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By solving the equilibrium equations for a polymer in a neutral polymer brush, the degree of 
interpenetration of two polymer brushes in contact and near contact is calculated. These results are 
used to calculate values of the force of static friction in agreement with recent friction measurements 
for polymer brush lubricated surfaces. It is shown that at sufficiently light loads polymer brush 
coated surfaces can slide, with the load supported entirely by osmotic pressure, at a sufficiently 
large spatial separation so as to avoid entanglement, and hence static friction. 

PACS numbers: 68.35.Af,81.40.Pq,46.55.+d 
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I. INTRODUCTION 

A polymer brush consists of a surface with a fairly con- 
centrated coating of polymer chains, each one of which 
has one end tightly bound to the surface. They serve 
as extremely effective lubricants, producing friction co- 
efficients as low as 0.001 or less! Polymer brushes are a 
promising way to reduce friction to extremely low values. 
They have the disadvantage, however, that they must be 
immersed in a liquid solvent in order to function as a 
lubricant. Since joints in the human body are known 
to be immersed in a fluid, known as the synovial fluid, 
it is likely that they are lubricated by polymer brushes, 
where polymers known as hyaloronan or lubricin are at- 
tached to the cartilage which coats the bones making 
up the jointsQ. The density profile of a polymer brush 
(i.e., the density of monomers as a function of distance 
from the surface to which the polymers are attached) is 
well cstabhshed0, What is not understood is how 
the interaction of polymer brush coated surfaces in con- 
tact with each other is able to account for the observed 
friction. For example, molecular dynamics studies gener- 
ally do not predict static friction f^, whereas surface force 
apparatus measurements due to Tadmor, et. al.j^, find 
that there is static friction. 

Pincus, Witten, Milner, etc. 0, 13 , present the fol- 
lowing picture: over all except the very outer edge of 
a brush, the monomer density at a distance z from the 

but there is a 
c is the height 



solid surface is proportional to h. 
low density tail for z > hmax, where hr, 
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of an uncompressed polymer brush. When two brushes 
are put in contact the brushes are compressed and the 
resulting elastic force (which is really of entropic origin) 
and the osmotic pressure support the load. The osmotic 
pressure is not of the usual textbook example in which a 
solute is trapped in an impenetrable membrane because 
here what plays the role of the membrane are the forces 
that hold the monomers together in each polymer chain. 
Since each polymer is anchored to a one of the two sur- 
faces, it cannot escape from the region between the two 
surfaces. Solvent is pulled into the region between the 
plates, diluting the monomer density (which plays the 
role of a solute density) resulting in osmotic pressure. 



The experiments of Tadmor, et. al. |^ exhibit static 
friction which increases logarithmically with the time 
that the polymer coated surfaces are in contact. The 
logarithmic increase of static friction with resting time 
is qualitatively similar to a similar effect for solid-solid 
contact, even though the mechanisms for it are much 
different 1^]. 

Kinetic friction is discussed in a paper by JoannyQ, 
assuming that it occurs by reptation of polymers from 
one brush that get entangled in the second brush. He 
finds only a viscous friction, but if the polymers near the 
interface get entangled enough, there must be something 
that resembles static friction if the times scale associ- 
ated with the reptation motion is very long compared to 
the time for the surfaces to move some very small dis- 
tance in creep motion as they attempt to slide. Even for 
solid-solid friction, there is always creep motion. Creep, 
however, for the polymer brushes happens over a much 
shorter time scale than for solid-solid contact. 

In order to calculate the static friction it is necessary 
to determine the number of monomers belonging to poly- 
mers of one brush that penetrate into the second brush. 
This can be accomplished using an extension of the ap- 
proach used in Ref. 0. Consider two polymer brushes 
which are pressed together. The mean field potential 
for the hard core interaction of a monomer belonging to 
a polymer under consideration with other monomers is 
given by 



w[(t>{z) + ^{D ~ z)], 



(1) 



where D is the separation between the surfaces, to which 
the two polymer brushes that we are considering are at- 
tached, z is the distance from one of the two surfaces, 
which we will consider to be the lower surface, (l){z) is the 
monomer density of one of the polymer brushes and w 
is the standard monomer excluded volume parameterjp . 
In equilibrium, has the standard parabolic density 2] 



w(l){z) = B{k 
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for z < h and zero for z > h, where h is the actual height 
of the brush and where hmax is its maximum height, 
which is the equilibrium height when the brush is not in 
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contact with a second brush and B = TT^ksT / {^a^N"^)^ 
the form that B must have in order for a polymer at- 
tached on one end to the surface at z=0 to extend to an 
arbitrary height z = p < h. In Ref. 0, it is shown that 
h^ax = (12/7r2)i/3(aVu;/fcBT)i/37V, where the factors 
of ksT and a, which were set equal to one in Ref. 0, 
have been inserted here. Here, a is the intermonomer 
distance, fc^ is Boltzmann's constant, T is the absolute 
temperature and a is the number of polymers belong- 
ing to the polymer brush attached to the surface per 
unit area. When the brushes are in contact and the two 
brushes are pushed together under a load, h < hmax, 
which means that the brushs get "flattened out" so that 
(j){z) no longer goes to zero at z=h. For infinite N, all 
of the monomers belonging to the polymers belonging to 
the lower brush (i.e., the polymers attached to the lower 
surface) remain within that brush (i.e., in the region for 
which z < h), and the monomers belonging to the poly- 
mers of the upper brush (i.e., the polymers attached to 
the upper surface) remain in that brush or in the region 
defined hy h < z < D, where D=2h). For finite values 
of N, however, there exist thermal fluctuations in which 
polymers from the lower brush penetrate into the upper 
brush and polymers belonging to the upper brush pene- 
trate into the lower one. In order to calculate the degree 
of penetration, let us consider the free energy needed to 
add a polymer to one of the brushes ^2], 

dr 

F{k)^F{k~l)= / dn[{l/2)(kBT/a^)\-^\'+ 
Jo dn 



w{cj){Zn) -t- cj){D - Zn))], 



(3) 



where F(k) is the free energy for a two brush system 
containing k polymers and r„ is the location of the n*'' 
monomer in the polymer that was added to the lower 
brushy]. We can find the a local minimum value of the 
single polymer free energy in Eq. (3) by setting the vari- 
ation of this expression with respect to Tn equal to zero, 
which gives 



{kBT/a^ 



dv? 



d(j){zn) d(j){D — z) 



dz„ 



dz„ 



(4) 



and dxn/dn — dyn/dn are constant. Since we are in- 
terested in small fluctuations from the infinite N con- 
figuration of the brushes (for which the brushes do not 
interpenetrate), to lowest order in 1/N, we may use the 
parabolic monomer density of Eq. (2), which was found 
in the infinite N limit. Solutions to this problem for which 
the polymer under consideration, which is attached to the 
lowest surface, penetrates into the region for which z > h 
may be found by finding such solutions to Eq. (4). For 
z < h, Eq. (4) is 



d(j){z) 



d^ 
dv? dz 
and for h < z < it \s 



d^ 
dn^ 



w- 



d(j){D 



dz 



-2Bz 



-2B{z-D). 



(5a) 



(56) 



The solution of Eq. (5a) is 

z = ccos{Ljn) (6a) 

and the solution to Eq (5b) is 

z = D + acos{Lon) -t- bsin{ujn) (66) 

where a, b, and c are constants and the "frequency" 
uj = [2Ba? /ksTf/'^ = ■n/{2N). The choice of the cosine 
solution in Eq. (6a) was made because for a polymer 
attached to the lower surface z^v = and the expression 
for UJ under Eq. (6b) will result in this condition being 
satisfied if we choose the cosine solution. If we assume 
that for n > ki, the value of n for the monomer which 
first penetrates into the second brush, we must demand 
that the values of z given by Eq.'s (6a) and (6b) and their 
derivatives with respect to n must be equal at n = fci, 
since these are second order differential equations. This 
condition gives 

a — ~h[cos{u;ki) — sin(ujki)tan[ujki)], (7a) 



-2hsin{ujki), 



h 

cos(ujki) 



(76) 



(7c) 



Since we expect for large N and low temperature that 
ki << N, we find from Eq. (6b), using Eq. (7) that 



{3/2)h{ujkif 



and 



^U=o ~ -2huj'^ki. 
dn 



(8a) 



(86) 



Following the discussion in Ref. 0] , we find that the fluc- 
tuation in the free energy of a polymer belonging to the 
lower brush when fci of its monomers penetrate into the 
upper brush is equal to the amount of work by the ten- 
sion that would have to be applied to the n=0 monomer 
to pull it out of the lower brush by an amount y — ZQ — h, 
which is given by 



fh+V 

AS = (fcsT/a^) / 

Jh 



dZn 

dn 



=odzo. 



(9) 



If we change variable of integration to fci using Eq.'s (8), 
we obtain 



AS* = 2{kBT/a^)Lo^h^kl 



(10) 



Using Eq. (10), we find that the probability of ki 
monomers of a polymer from the lower brush penetrating 
into the upper brush is given by the Boltzmann factor 



-AS/KbT 



exp[~2uj^{h/afkl^ 



(11) 
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Thus, we conclude that the mean length of a segment of 
a polymer belonging to one brush that gets entangled in 
the second brush is < /ci > found from Eq. (14), which is 
of order [{l/2){a/hraaxY'^~'^Y^^ oi' larger (for h < hmax, 
which is of order iV^/^, and hence these segments are 
sufficiently long for large N to form "blobs" which are 
close packed in the brush that it has penetrated. This is 
expected to occur because the polymers belonging to the 
two brushes are in the semi-dilute regime |^. This must 
be the case because the Flory radius i?_F of a polymer 
belonging to a polymer brushy is always greater than the 
mean spacing of the polymers, and hence the polymers 
belonging to the brushes cannot individually curl up into 
non-overlapping "blobs" of radius Rp, as they would in 
the dilute limit. The expression for the monomer density 
in a polymer brush ^(z) given by Eq. (2), can be written 
as (t){z) = iNa/hmax)[i ~ (z/hmax)'^]- The mean value of 
the monomer density in the region of the brush into which 
polymers belonging to the brush in which it is in contact 
have penetrated is very close to c = (j){z — h). Then, 
from Ref . , the mesh size in the brush is approximately 
i = a(ca3)-3/4 ^ 88a = 2.64 x IQ-^cm. Then the static 
friction force f per polymer that penetrates is the force 
needed to pull a "blob" through the mesh, which is 



/ 



1.52 X 10~*dyn, 



(12) 



which gives a force per unit area of J a ~ 2.15 x KfiN/m?. 

Let us now consider the exciting possibility that at suf- 
ficiently light load the bulk of the two polymer brushes 
might not be in contact, meaning that D > 2/i, and that 
the load is entirely supported by osmotic pressure due 
to those polymers that extend out of each brush into a 
thin interface region of thickness D — 2h. (Because of 
the osmotic pressure in the interface region, h will be 
compressed below h^ax-) K will be shown that this in- 
terface region can be sufficiently thick so that the poly- 
mers from the two brushes do not entangle to any sig- 
nificant degree, resulting in negligibly small static fric- 
tion. Then, the end of the polymer zq is located in 
the interface region and for the part of the polymer in 
the interface region (0 < n < ki), dzn/dn = — fi, a 
constant, and hence z„ = zq — vin. For N > fci, Zn 
is still given by Eq. (6a) and (7c). Again requiring 
continuity of Zn and dzn/dn at n = fci and assum- 
ing that Loki << 1, we obtain vi — hujki, and hence, 
AS* = [kBT/a^) j^' vidzQ = {2/i){kBT /a^)h'^uj^kl and 
hence the probability that fci monomers of the chain pen- 
etrate into the interface region is given by 



P ^ exp[~{2/2,){hLo'^ /afkl] 



(13) 



Hence < ki >w (3/2)i/3(a//iw^2/3 ^ \,12N'^/^ for the 
values for a and N from Ref.j^] used here and taking 
a ~ 3 X 10~^cm. In order for there to be no significant 
penetration of polymers belonging to one brush into the 
second, we require that D — 2h> zq — h ~ hxJ^ < ki >^. 
Substituting for < fci > from the expression under Eq. 



(13) for it, using the parameters from Ref. we obtain 
Zq — h ^ ATA". We can estimate the monomer den- 
sity in the interface region from c « [<^/{D — 2h)] < 
ki >« [(t/{zo - h)] < fci >« 0.4 X lO^^CTO^^. In the 
interface region for the parameters used here, we are in 
a dilute regime, in which the osmotic pressure is given 
by[H {l/2)kBTa^c^ for the above value of c is equal to 
0.864 X 10^A''/m^ at room temperature. Thus, at loads 
smaller than 0.864 x 10^7V/m^ the polymer brush coated 
surfaces would be separated by a sufficiently thick inter- 
face region that there would be negligible entanglement 
of polymers belonging to the two brushes, and hence 
there would be negligible static friction. Charged poly- 
mer brushes have counterions which can provide an addi- 
tional source of osmotic pressure, in addition to the net 
charge of the lubricated surfaces to support part of the 
load. This opens up the possibility that the load might be 
completely supported by osmotic pressure and electrical 
charge for charged polymer brushes. 

The normal force per polymer that can be supported 
by two brushes, each with a parabolic monomer distri- 
bution is the derivative of the free energy per polymer of 
one of the brushes given in Eq. (33) of Ref. 2], which 
can be written as 



iVVu; ,[0.9 + 1.5A3 



/?2 

rnax 



N 



(14) 



where A — {hmax ~ h)/hmax- The resulting normal force 
Jn per unit area, which is equal to minus the derivative 
of F with respect to h multiplied by cr, is given by 



(9/2)(-^)A2. 



/l2 

max 



(15) 



In the experiments of Ref. "5*1 , the surface coverage den- 
sity of the polymers in the brushes is (84A°)~'^ = 1.42 x 
lO^^cm^^. we determine that h^ax = 451 x lQ~^cm from 
which we find that /^cr = 0.811 x lO^A^iV/m^, which 
is of the same order of magnitude as the load in Ref.|^. 
Combining this and Eq. (12) and the discussion above it, 
it is found that the coefficient of static friction is greater 
than 0.0265, in agreement with Ref. 0. If the surfaces 
to which the polymers are attached are perfectly smooth 
(e.g., for mica), / oc ■ As discussed earlier, we ex- 
pect there to be creep for polymer brushes which occurs 
over a much shorter time scale than for typical lubricated 
surfaces. Let us consider what would happen if we ap- 
ply a shear force below the maximum static friction to 
the two lubricated solids. On the average, an entangled 
polymer will work itself free in a time r, the reptation re- 
laxation timeQ, resulting in the center of masses of the 
two surfaces sliding a distance Ax/Ne relative to each 
other, where Ax is of the order of the spacing between 
attached polymers in the brush and is the number 
of entangled polymers. Dividing this by r and multiply- 
ing by Ne, we get a mean creep velocity of Ax/{t). For 
the polymers of Ref. |^, each of which contains about 
1300 monomers, r is found0 to be about 0.02 s. As Ax 
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in Ref. Q is about 84A°, we obtain a creep velocity of 
4.2 X 10 ~^cto/s, so that the surfaces will slide by about 
1.5 mm in an hour. 

The kinetic friction in the slow speed sliding limit, 
which must be comparable to the static friction can be 
found from the Tomlinson model 7], as follows: We may 
crudely treat an entangled polymer as a spring, one end 
of which is moving in a potential possessing several min- 
ima, taken to represent the potential due to polymers in 
the brush in which it is entangled. The opposite end of 
the spring is assumed to move with a sliding velocity v. 
Then the total potential of this polymer is given by 

V = {l/2)a{x-vt) + Vix), (16) 

where a is the force constant of the spring, V{x) is a 
potential which possesses several minima (which repre- 
sents the energy needed to pull a polymer out of a brush 
in which it is entangled, and t is the time). If a < 

2^1^^ \max, this potential possesses multiple minima0- 

Here ^Jf'^ \ max is the maximum force constant of the 
bottom of the entanglement potential well. Following 
the treatment of a semi-dilute solution of polymers in 
Ref. 0, the entanglement potential represents the en- 
ergy that must be expended to pull a "blob" through 
the mesh created by the other polymers among which 
the polymer under consideration is entangled. In order 
to pull the polymer through the mesh, each "blob" must 
be stretched so that it will become small enough to fit 
through. It follows that the force constant of the mesh 
potential must be equal to the force constant for a single 



"blob." Since the chain as a whole consists of a bunch 
of " blobs" connected in series the effective force constant 
of the chain must be equal to the force constant of a 
single "blob" divided by the number of "blobs" of the 
polymer chain under consideration which are entangled 
in the other brush. Then clearly the above condition for 
multistability is always satisfied. As the chain is pulled 
through the potential a potential minimum containing 
the end of the chain periodically becomes unstable, al- 
lowing the end of the chain to drop into a potential well 
of lower energy. When this occurs, the resulting drop in 
potential energy gets converted into kinetic energy, which 
is assumed to get rapidly dissipated into excitations of the 
system. Setting the rate at which such energy dissipation 
occurs equal to the product of the kinetic friction and v, 
we determine that there must be kinetic friction in the 
low V limit. 

We conclude that entanglement of polymers belonging 
to one brush in the second brush is able to account for 
the static friction observed in Ref. |^ for two polymer 
brushes in contact. Although it was estimated that there 
is creep of about 1.5 mm per hour, on short time scales 
there does appear to be true static friction. Osmotic 
pressure alone was found to be able to support a load 
of the order of lO^Pa while keeping the polymer brushes 
sufficiently far apart to prevent the entanglement that 
results in static friction. 
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